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ABSTRACT 



With the increasing application of power electronics to shipboard electric power 
systems, there is a need for studying power systems in an object-oriented, graphical 
environment in order to examine the interrelations of the various system components. 

This thesis presents a method of analyzing shipboard electric power systems in a 
graphical environment The simulation method used allows for real-time observation of 
the various power system measurements of interest 

The machine equations, including those for voltage, torque, and rotor angle are 
developed which permit analyzing the machine in any chosen reference frame depending 
on the study being undertaken. Regulators for both speed and field voltage are developed 
and the stability of the speed regulator is examined. A means of modeling a load for a 
synchronous generator is then shown which accurately interrelates the loading of the 
generator and the frequency and voltage output of the machine. This load is then 
connected to the synchronous generator and two different scenarios are examined, 
including a fifty percent load change and a fault study. Graphs showing the output of 
these studies are presented. 
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I. 



SYNCHRONOUS MACHINE MODEL DEVELOPMENT 



The primary means of shipboard electric power generation is with synchronous 
generators driven by steam turbines. The electrical and electromechanical behavior of the 
synchronous generator can be accurately predicted from the equations for a three-phase 
salient-pole synchronous machine. 

The rotor of the synchronous generator is equipped with a field winding and one or 
more damper windings. The rotor windings usually all have different electrical 
characteristics. Besides the differing electrical characteristics, the salient-pole machine 
exhibits a magnetically asymmetrical geometry due to the varying air-gap path for the 
magnetic flux. This magnetic asymmetry creates time-varying self- and mutual-inductance 
terms in the voltage equations for the stator variables. By transforming the equations into 
a reference frame which rotates at the same angular velocity as the rotor, the magnetic flux 
path remains constant, thereby eliminating the time-varying inductance terms. 

This chapter develops the voltage and electromagnetic torque equations for the 
synchronous generator in terms of machine variables and then converts these equations to 
the rotor reference frame via Park’s transformation [1], The developed model assumes a 
magnetically linear system thereby neglecting the effects of saturation. Saturation could 
be accounted for with the addition of another simulation segment for the model, however 
saturation data for the machines studied was unavailable. 
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A. VOLTAGE EQUATIONS EM MACHINE VARIABLES 



Voltage equations for the synchronous generator will be developed from a simple 
two-pole, three-phase, wye-connected, salient-pole synchronous machine shown in Figure 



with sinusoidally distributed windings, displaced 120°, with N s equivalent turns and 
equivalent resistance r s (subscript s denoting stator variables). A field winding and three 
damper windings (fd, kd, kql, and kq2 respectively) are installed on the rotor. All 
windings are assumed sinusoidal with applicable N equivalent turns and resistance r 
(subscript denoting particular circuit). The quadrature axis ( q axis) and direct axis (d axis) 

are used with reference to the magnetic axis of the rotor circuits. Voltages of the rotor 
damper windings (v w , v^,, and v^ 2 ) are for monitoring purposes only, the windings 

themselves are physically shorted. 

Voltage equations in machine variables [2] may be expressed in matrix form as 



1. The number of poles of the particular machine can be easily changed as this is one of 
the variables of the resulting torque equations. The stator windings are assumed identical 




( 1 ) 

( 2 ) 



where 




dbcs s 




(3) 



2 



q axis 




Figure 1. 2-Pole, 3-Phase Synchronous Machine. 
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and 



L = 
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( 6 ) 



In (1) and (2), the time differential operator is p. In (4) and (6), the subscript l denotes a 
leakage inductance term. The subscripts skql, skq2, sfd, and skd denote mutual 
inductances between stator and rotor windings. 

The magnetizing inductances are defined as 



^=\(L a -L,) (7) 

iw =!<£,, + £,) (8) 



where the term incorporating L4 refers to the inductance term which is independent of 
rotor position and the term with Lg is the inductance term which is position dependent. 
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The rotor variables are then referred to the stator windings via the appropriate turns 
ratio giving the following equations 
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where the primed variables (i.e. ) refer the appropriate rotor variable to the stator 

windings. 

The voltage equations expressed in terms of machine variables referred to the stator 
windings are 



^ abcs 




r s + P L, 
!'«>' 



p K 

<+pK 



^ abcs 




( 12 ) 



Note that the voltage equations (12) are valid for positive direction of stator current 
assumed out of the stator terminals. If motor action is assumed, the sign preceding 
must be changed. 
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B. TORQUE EQUATION IN MACHINE VARIABLES 

The field energy stored in the coupling field of a synchronous machine can be 
expressed as 



W, = - L„I)i, te - (i. to ) T I4i;* 






(13) 



where 



L&- — diag^Lj^j 2 Ly d 1^] (14) 

The electromagnetic torque is found by taking the partial derivative of the coenergy with 
respect to the rotor angular displacement 0 r . The magnetic system is assumed linear, 
therefore the coenergy equals the field energy, and upon taking the partial derivative, the 
electrical torque is 

r, = ( Y , {-j(U )r ^[ L ,-^ I ]i A ,+ ( w )r g|-[ L ;]iu} (15) 



where P is the number of poles of the machine. The torque is positive for generator action 
with the positive direction of stator current assumed out of the stator terminals. 

The torque and rotor speed are related by 

T,=-J(,j)P<K r +T, (16) 
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where J is the inertia of the machine in kilogram meters 2 (kg-m 2 ) or Joule seconds 2 (J*s 2 ). 



The input torque 7) is positive for a torque applied to the shaft of the synchronous 
machine. 

C. VOLTAGE EQUATIONS IN THE ROTOR REFERENCE-FRAME 

The machine voltage equations (1) and (2) are seen to contain flux-linkage terms 
which are dependent on rotor position 0 r . The analysis of the machine behavior is greatly 
simplified by transforming these time-varying inductances to constant inductance terms. 
This can be realized if the stator variables are transformed to the rotor reference frame. 
R.H. Park was the first to incorporate a change of variables in the analysis of synchronous 
machines [1]. 

A change of variables which formulates a transformation of the 3-phase variables of 
the stationary circuit elements (thus the need to refer the rotor variables to the stator) to 
the rotor reference frame may be expressed as 




(17) 



where 




(18) 
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( 19 ) 



and 
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( 21 ) 



The term K' denotes the transformation matrix of stator variables (the subscript-5) to the 
rotor reference frame (the superscript-r). It can be shown that the voltage equations (1) 
and (2) are transformed into 



VqdOs ~ r jqdOs + ®r^dqs + P^qdOs 

_./r _/;/r . / r 

V qdr ~ r r l qdr + P^-qdr 



where 



Wdqs) 



-K 

o 



( 22 ) 

(23) 



(24) 
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The flux-linkages can be expressed in the rotor reference frame by transforming (9) to 
obtain 



V 

^ qdOs 
l ' r 

K qdr 



KK 



:(L;) r (K;r' l; 



-l 



qdOs 
tr 
l qdr 



(25) 



Expanding the voltage equations given by (22) and (23) gives 



Vr q s=- r / q s+^rK+PKs 


(26) 


V* ®rKs + P X ds 


(27) 


v Q s = ~ r Jos + PK 


(28) 


/r / */r , y n 
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(29) 
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V kd = r kd^kd P^kd 


(32) 


and the expanded flux linkage terms become 




= — As 1 ?* + Am? ( iqs + 4?1 A?2 ) 


(33) 
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II 
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(37) 
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(38) 
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(39) 
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Expressing the voltage and flux linkage equations in terms of reactances rather than 
inductances gives the voltage equations in matrix format (40), where the expression for 
the applied field voltage is introduced 



y 

/r _ /r 
e xfd ~ V fd , 

jd 



(41) 



and the leakage and magnetic reactances are combined as 



X q = X h + X mq 


(42) 


Xd — + Xnj 


(43) 


X'w-X^ + X" 


(44) 


X' kq2 = X' Ucq2 + X mq 


(45) 


X'fd = x; f d + x* 


(46) 


X' kd = X' lkd + X md 


(47) 



The zero sequence terms have been omitted since balanced conditions are assumed. 

The model being developed will be connected to a non-infinite bus. This implies that 
the voltage and frequency do indeed vary with load, based on the droop characteristics of 
the particular regulator. With this in mind, it is necessary to solve for currents based on the 
bus voltages, therefore the currents are selected as independent variables. Rewriting (40) 
with the current and current derivative terms separated. 
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(48) 

The state equation format of (48) is used for the simulation, so the current derivatives 
must be solved for. As an example 



v = [ A]i + [B]p i 
pi = [B] _1 {v - [A]i} 



(49) 



and, with this format, the current derivatives are solved as shown in (50) 
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Denoting the inverse matrix elements as [B]' 1 =[b] = 



^11 ^12 
^21 



'61 



the desired differential current equations are obtained 
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Pi'fd = M*4 + rf ds - co r i; +co r ^ i" + co r ^i^] 



co. 
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Note the similar form of the <?- and the d-axis integral current equations. These will 
greatly simplify the SIMULINK [3] model. 

D. TORQUE EQUATION IN SUBSTITUTE VARIABLES 

The torque equation (15), transformed into rotor reference frame variables is 



% = (f )(f + i’i + iZK - V-K. + + %)«,] 



(57) 



and writing in terms of reactance 
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gives 



z = (^)[^(-/; +/; +&)/; -^-h; +£, +/; r 2 )4] 



■l)A[i 

‘2 2 ok 



(59) 



The rotor speed is obtained from (16) as 



P ^r=(-—)(T e -T l ) 



(60) 



E. ROTOR ANGLE 

The rotor angle of a synchronous machine is a useful parameter to monitor the torque 
production of the machine. It relates the angular displacement of the rotor and stator 
magneto motive force axes. In equation form 

s = e, - e„ 

= f'[o>,(5) - -MUR + 9,(0) - 6„(0) <61) 

where 5 is the rotor angle (in radians), and £, is a dummy variable of integration. 

Equation (61) can be implemented by solving equation (61) in terms of the available 
parameters. The rotor position, 0 r is available by integrating the rotor speed. However, 
the voltage angle 0 CT must be solved for as follows 

e„ =fco„<5) 4+e„(0) (62) 
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and 



= v g* cos 0 r + v* sin 0 r 

= J v l + vi cos[0, + tan" 1 (-^-)] 

v ? * 

= -J v l + vl cos[ J (D r dt + 0 r (0) + tan" 1 (^-)] 
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(63) 



therefore, from (62) 



e„ = £«>„(?) 4+e„( 0 ) 

= f to, A + 8, (0) + tan' ' (— ) (64) 

J v 



giving 

5 = 0 r -©„ = tan” 1 (— ) (65) 

v„. 



where a positive rotor angle is assumed for generator action. 
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n. SPEED AND VOLTAGE REGULATORS 



For the machine model to accurately simulate a shipboard power system it is 
necessary to include controllers for both the field voltage and the prime mover torque. 
The two principal controllers of a turbine generator set are the automatic voltage regulator 
and the speed governor. The automatic voltage regulator model consists of voltage 
sensing equipment, comparators and amplifiers controlling a synchronous machine which 
can be generating or motoring. The speed governor may be considered to have similar 
equipment but in addition it is necessary to take the turbine into account. 

A. SPEED REGULATOR 
1. Model Development 

The block diagram of a steam turbine and associated governor for a typical 
naval ship's service turbine generator (SSTG) is shown in Figure 2 [4], 

The reference speed input is the no-load frequency selected by the operator of 
the electric plant. The gain K is used to adjust the speed droop of the regulator. x s is the 
time constant associated with the servo and steam valve system andx p is the time constant 
associated with the conversion of steam to torque (T pm ). The following values were used 

for this study 

0.1 sec 

x p 0.06 sec 

K 150 

From this block diagram (Figure 2), a SIMULINK model was developed and is shown in 
Figure 3. 
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Figure 2. Block Diagram, Steam Turbine & Governor. 
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Figure 3. SIMULINK Speed Regulator. 
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As can be seen, the speed regulator is designed to compare the reference speed with the 
speed being produced by the synchronous machine. This speed difference determines the 
prime mover torque being applied to the generator. When the speed difference is no 
longer changing, the prime mover torque is equal to the electrical torque developed by the 
generator (load dependent). The final machine speed (and therefore the machine output 
frequency) will be at some value below the no-load set-point frequency. This can be 
demonstrated by referring to the 'house' diagram shown in Figure 4. 

As the output power of the machine is increased, the operating frequency of the 
machine decreases (based on the slope or droop of the particular regulator) the operating 
frequency can be restored to a desired frequency (i.e., 60 Hz) by the operator adjusting the 
no-load setpoint of the regulator. 

2. Regulator Stability 

As was mentioned in the previous section, the speed droop characteristics of the 
regulator are determined by the gain K of the regulator. The percent speed regulation is 
defined by 



co ,-co„ 

% Speed Regulation = — 100% (62) 

As the gain K is raised, the percent speed regulation approaches 100%, at which 
point the regulator and therefore the regulated plant will become unstable. This stability 
point was determined by examining the root locus of the model of Figure 3. The 
synchronous machine model was obtained by linearizing the machines current block 
around its stable operating point for rated conditions by using the MATLAB [5] function 
linmod. The linearized plant model returned by linmod is in linearized state space format 
which is then changed to transfer function format with the function ss2tf. The transfer 



21 




<U -g 
CO © 
cl, 



Figure 4. House Diagram. 
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Operating 

Power 



function format of the plant, speed block, and speed regulator is then reduced to a single 
open loop to study the roots of the system as the regulator gain K is varied. This 
complete process is demonstrated in Figures 5, 5a, 5b, and 5c. The resulting root locus 
plot is shown in Figure 6. The operating constraint was determined for the first root 
meeting the imaginary axis for a given K. The limit of K was determined to be 153. 
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Figure 5. Synchronous Machine Model. 
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Figure 5a. Current Block. 
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BLOCK VI 





Figure 5b. Linearized Current Block, Speed Regulator, and Speed Block. 
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Figure 5c. Transfer Function Development, Speed Regulator. 
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Figure 6. Root Locus Plot-Synchronous Machine and Speed Regulator. 
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B. VOLTAGE REGULATOR 



The voltage regulator or excitation system used in this study is shown in block 
diagram form in Figure 7 [6]. The parameters are defined as follows 

x A and K a Regulator amplifier gain and time constant 

x E and K e Saturating exciter gain and time constant 

x n , x p 2 , and K f Feedback loop gain and time constants 

The SIMULINK model developed from the block diagram is shown in Figure 8. The 
regulator compares the bus voltage with the specified reference voltage. As can be seen in 
the Simulink model, this comparison is done in the qdO reference frame. 

The reference voltage is computed as follows 

Vreference^^ + V/ (63) 



This equates to the desired line-to-neutral voltage (peak) as shown below 






(64) 



V 


2 


COS0 r cos(0 r ~~~) COS(0 r + ~) 
^ 2k n 2 k 


V 




~ 3 


sin0 r sm(0 r — — ) sm(0 r +— ) 




- v 0_ 




1 1 1 


_ v c. 






.2 2 2 





(65) 
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Figure 7. Block Diagram, Voltage Regulator. 
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Figure 8. SIMULINK Voltage Regulator. 
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F1 V F2 



where 



v a =v M cos(0 J 
v b =v i n cos(0 ev 

2k 

v c =v in cos(e £V +— ) 



and 



e r =J o 'm r ©y+e,(0) 

e„=|>,©^+e„(0) 

and solving for Vq and 



( 66 ) 



(67) 



2 

V = — V; 

9 3 h 



COS0 r COS0 CT +COS(0 r — ^)COS(0 CT — ~~)+COS(0 r +~~)COS(0 ev +-J-) j 



( 68 ) 



2 

1 

3 



v,i=-v h 



sin 0 r cos0 CT +sin(0 r -y-)cos(0 ev -y)+sin(0 r +^y)cos(0 ev +-y)j 



(69) 



and using the trigonometric identities 



, 271. . 271. . 27t. . 2n. 3 , . 

cos j: cos y+cos(x )cos(y )+cos(xH )cos()H )=— cos(x-y) 

3 3 3 3 2 

. . 2ti. . 27t. . , 2tt. . 27t. 3 . . . 

sin xcosy+sin ( x )cos(y )+sm (xH )cos(y+ — )=— sin (x-y) 

3 3 3 3 3 



(70) 

(71) 
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gives 



V? = 3 Vfa 2 COS ( 0r ~ 0ev) 


(72) 


= v M cos[(co r -a>,)r+0 r (O)-0„(O)] 




v,=|v.|sin(e r - e J 


(73) 


=v te sin[(co r -a),>+e,(O)-0 ev (O)] 




At synchronous speed co r =co e , and from (61) 5=0 r -0 CT , therefore 




v„=v M cos(S) 


(74) 


v d =v M sin(8) 


(75) 


and 

y reference = V V m[COS 2 (8) + Sin 2 (8)] = V M 


(76) 



33 



ID. MACHINE LOAD 



To accurately study the machine's performance it is necessary to be able to load the 
given machine in a realistic manner. The inductive nature of transmission lines and 
induction motors dictate a need for a resistive-inductive (RL) loading of the machine. 

This section develops the RL-load used for the synchronous generator and then the 
SIMULINK implementation of this loading. 

A. PASSIVE RL-LOAD 

The load used for this study was comprised of a passive resistive-inductive load. 
This loading is typical of ac motors, transmission lines, electric heating, and lighting. The 
load was modeled in afcc-variables in order to eliminate the need for additional 
transformations (abc -» qdO). The simulation uses voltage as the input and currents as the 
outputs. This requires the phase voltages to be made up from the network equations. The 
loads also use voltage as the inputs and current as the outputs. Based on this topology it 
is advantageous to assume a Y-connected set of large resistors are connected across the 
terminals of the synchronous generator [7] and the phase voltages at the terminals of the 
generator are computed from 



v„ = - R(i aG +i al) 


(85) 


V b = ~ RdbG ) 


✓—N 

oo 

On 


*c= ~ WcG+D 


(87) 



where i aG , i hG , and i cG are the phase currents flowing into the generator from the load and 
i bL , and i cL are the phase currents flowing into the load from the generator. These 
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currents are computed first as i q , i d , and i 0 and then transformed to i a , i b , and i c . This 
is represented in Figure 9. 

This load modeling is quite accurate as the resistance value R is made large (20-50 
pu). As the resistance is made large, the current through this branch becomes negligible. 

B. SIMULINK LOAD MODEL 

The passive i?L-load was modeled in SIMULINK by starting with the voltage 
equations (only developed for phase a, others similar) from Figure 10. 



v a= e a + L s pi a +r s i a 


(88) 


v a =r L +L L pi L 


(89) 


V a =R(i a ~i L ) 


(90) 



The differential equations can be implemented as shown in Figure 11. 

As can be seen, the terminal voltage v a is dependent on the output current i a . This 
voltage dependency is the difference between a finite and an infinite bus model of electric 
machines. An infinite bus analysis assumes no dependency between v a and i a . 

An attempt was made to make the load purely resistive. SIMULINK is unable to 
deal with this type of load however since this would incorporate direct feed-through 
between input and output of the load. This forms an algebraic loop and cannot be 
analyzed. The inductance of the transmission line overcomes this limitation. The purely 
resisitve load, connected by a transmission line (with inherent inductance) to the generator 
bus adds sufficient inductance to warrant modeling. This inductance then incorporates a 
time-delay between input and output due to the integral form of the load equations. 



35 




Figure 9. Generator Terminal Voltage Model. 
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Figure 10. Passive RL-Load Model. 
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Figure 11. Differential Equation Implementation of Voltage Equations. 
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rv. MACHINE MODEL PERFORMANCE 



The synchronous machines performance is examined in this section. The machines 
parameters are given and the machines response to both a load change and a fault 
condition are studied. 

A. MACHINE PARAMETERS 

The machine being simulated is a 4-pole, 60 Hz, 440 V (line-to-line), 500 kW, 0.85 
power factor synchronous generator. The per unit parameters on the machine base (Z 5 ^ 
= 0.1291 O) are as follows 

stator resistance, r s = 0.00031 

rotor field resistance, /yy = 0.00226 

q-axis damper winding resistance, 77 ^ = 0.0224 

d-axis damper winding resistance, rj^ - 0.00673 

stator leakage reactance, = 0.0632 

d-axis magnetizing reactance, X^ = 1 .4648 

q-axis magnetizing reactance, = 0.7668 

field leakage reactance, X/yy = 0.1545 

q-axis damper winding leakage reactance, X/£g = 0.077 1 

d-axis damper winding leakage reactance, X//^ - 0.3676 

rotor inertia, J = 14.526 

These machine parameters are taken from the Machinery Systems Laboratory (MSL) test 
facility for a typical shipboard synchronous generator [ 8 ]. The MSL test facility is within 
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the Machinery Systems Engineering Branch at David Taylor Naval Ship Research and 
Development Center. 

B. SIMULATION: LOAD CHANGE 

The machine was initially loaded at an approximate power factor of 0.8, with the load 
calculated as follows 

Machine data: 440-V 1.5 MVA . From Figure 12 



/I = L5e6 1 = 1. 969M 


(91) 


1 3 254 


/ = 1. 969 £--' eos ' I(0 - 8) m 


(92) 


Z = --Jlo - = 129 e ;36 ' 87 °mQ 

1.969e- j36g7 kA 


(93) 


r L = 103. 2mO X L = 77. 42 mQ. 


(94) 



With the machine loaded as above, the expected steady-state power delivered to the 
load will be approximately 1.2 MW. This does not however account for the losses in the 
stator windings. The machine then underwent a 50 % step load change by doubling the 
values of and . This was done in SIMULENK by the load-change block shown in 
Figure 13. The plots of the response are shown in Figures 14 through 21. The data 
obtained from the simulation was as follows: 





FULLY LOADED 


50% LOADED 




MODEL 


MODEL 


Load power, MW 


1.1754 


0.5186 


Machine power, MW 


1.1792 


0.5151 


Frequency, Hz 


59.99 


64.32 


Power Factor 


0.7992 


0.7832 
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Figure 12. Initial Load Parameter Determination Circuit 
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Figure 13. SEV1ULINK Load Change Block. 
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Figure 14. 50 Percent Load Decrease, vq. 
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Figure 15. 50 Percent Load Decrease, vd. 
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Figure 16. 50 Percent Load Decrease, va. 
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Figure 17. 50 Percent Load Decrease, field voltage. 
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Figure 18. 50 Percent Load Decrease, ia. 
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Figure 19. 50 Percent Load Decrease, te. 
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Figure 20. 50 Percent Load Decrease, rotor angle (del). 
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Figure 21. 50 Percent Load Decrease, rotor speed (co r/ ). 
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C. SIMULATION: THREE-PHASE FAULT 



The stability of a synchronous machine throughout a power system following a fault 
is of major concern. A 3-phase fault or short-circuit rarely occurs and a 3-phase fault at 
the machine terminals is even more uncommon; however, it is still instructive to 
demonstrate the ability to insert a fault for a chosen time interval and study the machines 
dynamic performance. 

The machine was initially loaded at a power factor of 0.8 as in part B. A 3-phase 
fault was inserted and was left in place for 7 cycles. The fault impedance (representative 
of the bus work) is 0.001291 + j 0.1291 Q. If the fault had been allowed to remain on the 
system slightly longer the machine would have become unstable in terms of not having 
returned to synchronous speed after removal of the fault or the machine would have 
slipped poles prior to returning to synchronous speed. It was determined that a seven- 
cycle fault is representative of a typical breaker reclosure timing cycle which the machine 
could have been subjected to. 

The three-phase fault was performed in SIMULINK by the fault block shown in 
Figure 22. The fault implemented by using a step function block which stepped in the 
fault impedance at a particular time and then another step function block stepped in the 
normal load at the end of the desired fault time. The plots of the response are shown in 
Figures 23 through 30. 
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Figure 23. 7-Cycle, 3-Phase Fault, vq. 
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Figure 24. 7-Cycle, 3-Phase Fault, vd. 
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Figure 25. 7-Cycle, 3-Phase Fault, va. 
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Figure 26. 7-Cycle, 3-Phase Fault, field voltage. 
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Figure 27. 7-Cycle, 3-Phase Fault, ia. 
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Figure 28. 7-Cyde, 3-Phase Fault, te. 
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Figure 29. 7-Cyde, 3-Phase Fault, rotor angle (del). 
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Figure 30. 7-Cycle, 3-Phase Fault, rotor speed ((D,.). 
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V. 



CONCLUSIONS 



A. ADVANTAGES 

The advantages to studying naval power system components in a graphical 
environment are many. To be able to see, in real time, any line or load parameter (i.e. 
current or voltage) makes the study of system component interaction very straightforward. 

The concept of "closing the loop" to make the input and output of the machine 
dependent on one another, as opposed to the infinite bus concept, is easily implemented 
with the graphical tools available. Navy power systems definitely require this non-infinite 
bus modeling due to their capacity difference as compared with commercial power grids. 

B. ANALYSIS OF SIMULINK 

Matlab's SIMULINK is an excellent tool for the type of studies conducted in this 
thesis. Model building and integration into the "power system" was easily done. 

One of the difficulties with SIMULINK included the algebraic loop problem 
necessitating an integrator in the path to provide a time delay element In this particular 
study this was not a problem due to the resulting differential equations, however, a 
possible purely resistive load would be troublesome. 

On the issue of speed of the simulations, naturally one is always looking for a faster 
simulation. I began my research on a 486-33MHz machine, then graduated to a SPARC-2 
machine and realized a speed difference of about two times. I then moved to a SPARC-10 
and achieved a speed increase of approximately four times over the SPARC-2. Typical 
times to run a simulation of approximately 5-seconds was about 4-hours on the SPARC- 
10. In comparing this to a non-graphical simulation tool such as Advanced Continuous 
Simulation Language (ACSL) [9], a comparable model runs approximately in real time (60 
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seconds of run time data in approximately 60 seconds). The disadvantage of ACSL is, of 
course, the inability of the user to observe system output in real time and the non-modular 
modeling of systems. 

C. FUTURE WORK 

The modeling of additional load components including induction machines, another 
synchronous generator sharing the load, parallel RL-loading, and adjustable speed drives 
could all be performed in SIMULINK. As the number of integrators in the complete 
model increases, however, the run time will also increase. In addition, the existing RL- 
load can be examined under such conditions as critical clearing time, first swing transient 
stability limit, and various fault conditions such as single line to ground faults, line-to-line 
faults, etc. 

In addition, voltage regulator stability can be studied to find gain limits for this 
particular regulator. Different regulator designs could also be examined. 

Mathworks has also advertised a software package which will take existing 
SIMULINK models and compile them; this will supposedly increase simulation speed 
considerably. This product sounds very encouraging. 
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